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ABSTRACT 


The  technique  of  image  compression  using  Iterative  Function  System  (IFS)  is 
known  as  fractal  image  compression.  An  extension  of  IFS  theory  is  Partitioned 
or  local  Iterative  Function  System  (PIFS)  for  coding  the  gray  level  images.  Several 
techniques  of  PIFS  based  image  compression  has  already  been  proposed  by  many 
researchers.  The  theory  of  PIFS  appears  to  be  different  from  the  theory  of  IFS  in 
the  sense  of  application  domain.  In  the  present  article  we  have  proposed  a  math¬ 
ematical  formulation  for  the  existence  of  the  attractor  of  PIFS,  assuming  it  as  a 
separate  scheme,  in  the  context  of  image  compression  .  It  has  been  shown  that  the 
attractor  is  an  approximant  of  the  given  target  image.  The  experimental  results 
have  also  been  presented  in  support  of  the  theory.  The  experimental  results  have 
been  obtained  by  using  a  GA  based  PIFS  technique  proposed  by  Mitra  et  al  [1]. 


1  Introduction 

The  theory  of  fractal  based  image  compression  using  Iterative  Function  System 
(IFS)  was  proposed  by  Barnsley  [2,  3].  He  modeled  real  life  images  by  means 
of  deterministic  fractal  objects  i.e.,  by  the  attractor  evolved  through  iterations 
of  a  set  of  contractive  affine  transformations.  Once  the  set  of  contractive  affine 
transformations  (say)  is  obtained  the  rest  is  an  iterative  sequence  of  the  form 
{J^^{0)}n>o^  where  ‘‘O”  is  an  initial  object  to  start  the  iterative  sequence.  The 
set  of  contractive  affine  transformations  is  called  IFS.  In  particular  at  the  Nth 
iteration,  the  object  On  is  used  as  input  to  the  IFS,  where  On  is  the  output  object 
obtained  from  the  (A^  —  l)th  iteration.  The  detailed  mathematical  description  of 
the  IFS  theory  and  other  relevant  results  are  available  in  [2,  3,  4,  5,  6]. 

Image  compression  using  IFS  can  be  looked  upon  as  an  inverse  problem  of  iterative 
transformation  theory  [7],  The  basic  problem  here  is  to  find  appropriate  contrac¬ 
tive  affine  transformations  whose  attractor  is  an  approximation  of  the  given  image. 
Thus  for  the  purpose  of  image  compression  it  is  enough  to  store  the  relevant  pa¬ 
rameters  of  the  said  transformations  instead  of  the  whole  image.  This  technique 
reduces  the  memory  requirement  to  a  great  extent.  But  the  question  is  how  to 
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construct  the  ti*ansformations  for  a  given  image.  A  fully  automated  fractal  based 
image  compression  technique  of  digital  monochrome  image  was  first  proposed  by 
Jacquin  [7,  8,  9].  This  technique  is  known  as  partitioned  [10]  or  local  [3]  iterative 
function  system.  The  partitioned/local  IFS  is  an  IFS  where  the  domain  of  appli¬ 
cation  of  the  contractive  affine  transformations  is  restricted  to  the  small  portions 
of  the  image  (subimages)  instead  of  the  whole  image  as  in  the  case  of  IFS.  In  PIFS 
the  given  image  is  first  partitioned  into  non  overlapping  square  blocks.  In  encod¬ 
ing  process,  separate  transformations  for  each  square  blocks  are  then  found  out  on 
the  basis  of  their  similarity  with  other  square  blocks.  These  blocks  are  large  in 
size  and  are  located  any  where  within  the  image  support  of  the  given  image.  As 
the  transformations  are  applied  partition  wise  the  scheme  is  known  as  partioned 
IFS.  Again  as  the  similarity  of  square  blocks  are  found  out  locally  it  is  also  known 
as  local  IFS.  Different  schemes,  using  PIFS,  have  been  proposed  by  several  other 
researchers  [10,  11,  1,  12]. 

The  theory  of  partitioned/local  IFS  appears  to  be  different  from  the  theory  of 
IFS  in  the  sense  of  restriction  of  the  application  domain  for  the  contractive  affine 
transformations.  So,  the  questions  are,  how  does  PIFS  produce  an  attractor  and 
how  does  it  become  a  close  approximant  of  the  given  target  image.  Generally,  PIFS 
is  considered  as  a  simple  extension  of  IFS  and.  it  is  assumed  that  the  Theoretical 
foundation  of  PIFS  is  same  as  that  of  IFS.  The  present  provides  a  mathematical 
formulation  for  the  existence  of  an  attractor  of  partitioned  IFS,  assuming  it  as  a 
separate  scheme,  in  the  context  of  image  compression.  In  particular,  firstly  we 
have  shown  that  the  transformations,  in  the  PIFS  scheme,  give  rise  to  a  fixed 
point  (attractor).  Secondly,  it  has  been  shown  that  the  transformations,  though 
not  exactly  contractive,  are  eventually  contractive.  Finally,  we  have  proved  the 
attractor  and  the  gi  ven  image  are  very  close  to  each  other  in  the  sense  of  a  chosen 
metric  or  distance  measure. 

In  the  next  section  we  have  described  the  theory  of  image  coding  using  IFS.  Sec¬ 
tion  3  consists  of  basic  features  of  constructing  PIFS  codes  for  a  given  image. 
Section  4  deals  with  the  basic  difference  between  IFS  and  PIFS  techniques.  The 
proposed  mathematical  formulation  of  PIFS  hets  been  discussed  in  Section  5.  The 
experimental  results  have  been  presented  in  Section  6  and  the  conclusions  are 
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drawn  in  Section  7. 


2  Theoretical  Foundation  of  Image  Coding  By 
IFS 

The  salient  features  of  IFS  theory  and  image  coding  through  IFS  are  given  below. 

Let  {X,  d)  be  a  metric  space,  where  X  is  a  set  and  d  is  a  metric.  Generally,  A  is 
taken  as  the  collection  of  compact  sets  and  d  is  taken  as  distance  measure  between 
two  sets  in  X.  Let  /  be  a  contractive  affine  map  defined  on  metric  space  {X,d) 
such  that  f  :  X  X  and  d(/(xi),/(x2))  <  ^  d{xi,X2)\  Vxi,X2  €  X,  where 
0  <  s  <  1  is  called  contractivity  factor  of  the  map  /.  For  any  large  positive 
number  N,  lim  f^{x)  =  a,  Vx  e  A,  and  also  /(a)  =  a.  “a”  is  called  fixed  point 

JV-'+oo 

(attractor)  of  /.  Here  f^{x)  is  defined  as 

f^{x)  =  /(  ),  with  f^{x)  =  f{x),  yx  e  X, 

Now,  let  /  be  a  given  image  which  belongs  to  the  set  X,  Our  intention  is  to  find 
a  set  T  of  affine  contractive  maps  for  which  the  given  image  /  is  an  approximate 
fixed  point.  T  is  constructed  in  such  a  way  that  the  distance  between  the  given 
image  and  the  fixed  point  (attractor)  of  T  is  very  small.  To  any  set  S  belongs  to 
A^,  the  set  of  transformations  T  is  used  as  follows; 

^(5)  =  U/i(5)- 

t 

The  attractor  “>1”  of  the  set  of  maps  T  is  defined  as  follows  : 

lim  T^U)  =  A,  VJ  G  A, 

iV-^OO 

and  where  ^^{J)  is  defined  as 

),with 

J^\J)  =  vjeA. 

Also  the  set  of  maps  J-  is  defined  as  follows; 

d{HJ^),HJ2))  <  s  d(Ji,  J2);  VJi,  J2  6  A  and  0  <  s  <  1.  (1) 
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s  is  called  the  contractivity  factor  of 


Let  d{I,:F{I))<t  (2) 

where  c  is  a  small  positive  quantity.  Now,  by  Collage  theorem  [2],  it  can  be  shown 
that 

d(I,A)  <  ~  (3) 

i  —  5 

where  A  is  the  attractor  of  T. 

From  (3)  it  is  clear  that,  after  a  sufficiently  large  number(iV)  of  iterations,  the  set 
of  affine  contractive  maps  T  produces  a  set  which  belongs  to  X  and  is  very  close 
to  the  given  original  image  7.  Here,  (X,  T)  is  called  iterative  function  system  and 
T  is  called  the  set  of  fractal  codes  for  the  given  image  I. 

In  the  context  of  digital  monochrome  image,  the  coding  scheme  is  called  partitioned 
or  local  Iterative  Function  System.  In  the  next  section,  the  construction  of  PIFS 
codes  has  been  described. 


3  Construction  of  PIFS  Codes 

The  structure  of  PIFS  codes  are  almost  same  as  that  of  IFS  codes.  The  only 
difference  is  that  PIFS  codes  are  obtained  and  applied  to  a  particular  portion  of 
the  image  instead  of  the  whole  image. 

Let,  /  be  a  given  image  having  size  w  x  w  and  the  range  of  gray  level  values  be 
[0,(?].  Thus  the  given  image  /  is  a  subset  of  IR^.  The  image  is  partitioned  into  n 
non  overlapping  squares  of  size,  say  6x6,  and  let  this  partition  be  represented  by 
TZ  =  {7^i,7?.2,  ’  •  •  Each  TZi  is  named  as  range  block.  Note  that  n  =  f  x  f. 

Let  V  be  the  collection  of  all  possible  blocks  which  is  of  size  26  x  26  and  all  are 
within  the  image  support.  Let  V  =  {Vi,T>2,  •  •  *  Each  Vj  is  named  as 

domain  block  with  m  =  (itr  —  26)  X  (u;  —  26)  and  Vj  C  [0,  u;]  X  [0,  u;]. 

Let, 

=  {/  :  X>j  — >  Dl^  ;  /  is  an  affine  contractive  map}. 

Now,  for  a  given  range  block  72.,*,  let,  fiy  G  JFj  be  such  that 

fium))  <  diiZi,  f{T>j))  V/  e  Vj. 
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Now  let  k  be  such  that 


diTli,  fi\k{Vk))  =  mm{  fnj{V^))  }  (4) 

Also,  let 

Our  aim  is  to  find  fi\k{T>k)  for  each  i  G  {1, 2,  •  •  • ,  n}.  In  other  words,  for  every  range 
block  TZi,  one  needs  to  find  an  appropriately  matched  domain  block  Vk  as  well  as 
an  appropriate  transformation  fi\k-  The  set  of  maps  ^  =  {/i|»)/2l»)  ■ ' '  ^/n]*}  thus 
obtained  is  called  the  partitioned  or  local  IFS  or  fractal  codes  of  image  7. 

To  find  the  best  matched  domain  block  as  well  as  the  best  matched  transformation, 
we  are  to  search  all  possible  domain  blocks  as  well  as  all  possible  transformations 
with  the  help  of  equation  (4).  The  Problem  of  searching  appropriately  matched 
domain  block  and  transformation  for  a  range  block  has  already  been  solved  by 
enumerative  search  [8]  and  by  using  Genetic  Algorithms  [1]. 

The  affine  contractive  transformation  /i|,  is  constructed  using  the  fact  that  the 
gray  values  of  the  range  block  are  scaled  and  shifted  version  of  the  gray  values 
of  domain  block.  The  contractive  affine  transformation  fiy  defined  on  IR^  is  such 
that  fi\j{Vj)  ->  Tli.  Also  fi\j  consists  of  two  parts,  one  for  spatial  information  and 
the  other  for  information  of  gray  values.  The  first  part  indicates  which  pixel  of  the 
range  block  corresponds  to  which  pixel  of  domain  block.  The  second  part  is  to  find 
the  scaling  and  shift  parameters  for  the  set  of  pixels  of  the  domain  blocks  to  the 
range  blocks. 

The  first  part  is  shuffling  the  pixel  values  of  the  domain  block  and  can  be  achieved 
by  using  any  one  of  the  eight  possible  transformations  (isometry)  on  the  domain 
blocks[8,  11].  Once  the  first  part  is  obtained,  second  part  is  estimation  of  a  set 
of  values  (gray  values)  of  range  blocks  from  the  set  of  values  of  the  transformed 
domain  blocks.  These  estimates  can  be  obtained  by  using  the  least  square  analysis 
of  two  sets  of  values  [11,  1). 

The  second  part  is  obtained  using  least  square  analysis  of  two  sets  of  gray  values 
once  the  first  part  is  fixed.  Moreover  the  size  of  the  domain  blocks  is  double  that  of 
the  range  blocks.  But,  the  least  square  (straight  line  fitting)  needs  point  to  point 
correspondence.  To  overcome  this,  one  has  to  construct  contracted  domain  blocks 
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such  that  the  number  of  pixels  in  the  contracted  domain  blocks  become  equal  to 
that  of  range  blocks.  The  contracted  domain  blocks  are  obtained  by  adopting 
any  one  of  the  two  techniques.  In  the  first  technique  the  average  values  of  four 
neighboring  pixel  values  in  a  domain  blocks  are  considered  as  the  pixel  values  of 
the  contracted  domain  blocks  [8].  In  the  other  scheme,  contracted  domain  blocks 
are  constructed  by  taking  pixel  values  from  every  alternative  rows  and  columns  of 
the  domain  blocks  [11,  1]. 

Now  to  select  appropriately  matched  domain  block  (X>jt)  and  appropriately  matched 
transformation  {fi\k)  for  a  range  block  (7?.,),  the  distance  measure  ‘‘d”  plays  an  im¬ 
portant  role.  The  distance  measure  ‘‘d”  [used  in  equation  (4)]  is  taken  to  be  the 
simple  Root  Mean  Square  Error  (RMSE)  [or  Mean  Square  Error  (MSE)]  between 
the  original  set  of  gray  values  and  the  obtained  set  of  gray  values  of  the  concerned 
range  block.  Let  Tliip^q)  ^i\k{p^q)  be  respectively  the  original  and  the  ob¬ 
tained  values  of  the  {p.qY^  pixel  of  the  range  block  Tli.  Thus,  the  expression  for 
RMSE  will  be 

d{1Zi,  Tluk)  = 

The  RMSE  is  not  a  metric  though  it  serves  the  purpose  of  a  distance  measure.  As 
selection  of  PIES  code  for  a  range  block  is  dependent  only  on  the  estimation  of  pixel 
values  of  that  block,  it  is  enough  to  calculate  only  the  distortion  of  the  original  and 
estimated  pixel  values  of  the  block.  The  first  part  of  the  PIES  code  indicates  pixel 
to  pixel  correspondence  of  the  contracted  domain  block  and  the  concerned  range 
block.  So,  this  part  of  the  transformation  is  fixed  for  every  iteration.  Only  the 
gray  values  of  the  range  block  are  being  changed  through  iteration.  Thus,  RMSE 
is  taken  as  the  distance  measure.  Note  that  the  same  measure  had  been  used  in 
all  most  all  the  articles  [8,  10,  11,  1,  12]. 

4  How  PIFS  technique  differs  from  IFS 

An  extension  of  the  iterative  function  system  concept  is  the  partitioned  iterative 
function  system.  PIES  mainly  differs  from  lES  in  the  domain  of  application  of 
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their  respective  transformations.  In  PIFS  the  transformations  are  not  applied  to 
the  whole  image,  as  in  the  case  of  IFS,  but  rather  have  restricted  domains.  In  all 
PIFS,  a  transformation  /,•  is  specified  not  only  by  an  affine  map,  but  also  by  the 
domain  to  which  /j  is  applied. 

The  difference  in  the  domain  of  application  of  the  two  techniques  is  shown  in  Figure 
1(a)  and  Figure  1(b).  Three  affine  contractive  transformations  are  applied  on  the 
image  Iq  to  result  in  an  image  which  consists  of  three  parts  /qi,  Iq2  and  I03  (Fig. 
1(a)).  These  three  transformations  are  then  applied  sequentially  to  result  in  a  fixed 
point  (attractor).  The  set  of  transformations  {/i,/2,/3}  is  called  the  IFS,  and  the 
attractor  of  this  set  of  transformations  in  this  case  is  Sierpinski  gasket  [2]. 

Contrary  to  the  above.  Fig.  ^(b)  is  showing  that  in  PIFS  the  map  /,•,,•  is  applied  to 
the  domain  Vj  result  in  which  is  an  estimate  of  71,-.  In  the  next  iteration, 
this  estimate  (71.)  is  not  used  as  the  input  to  the  map  fiy.  In  particular  in  the 
^xt  iteration  an  estimate  of  Vj  is  used  as  the  input  to  obtain  improved  estimate, 
T?.,-,  of  TZi.  A  domain  block  includes  many  other  range  blocks  or  part  of  them  (Fig. 

1(b)).  So,  the  estimate  of  Vj  consists  of  several  other  estimated  range  blocks  or 
part  of  them 

Another  important  and  significant  difference  of  PIFS  and  IFS  lies  in  the  context 
of  contractivity  factor  of  the  transformations.  For  an  IFS  with  an  expansive  map 
/,-,  the  set  of  maps  will  not  converge  to  a  compact  fixed  point.  The  expansive 
part  will  cause  the  limit  set  to  be  infinitely  stretched  along  some  direction.  This 
is  not  necessarily  true  for  a  PIFS.  PIFS  can  contain  expansive  transformations 
and  still  have  a  bounded  attractor.  So,  it  is  not  necessary,  in  PIFS,  to  impose 
any  contractivity  condition  on  all  the  transformations.  A  sufficient  contractivity 
requirement  is  that  the  set  of  transformations  T  be  eventually  contractive  [10]. 
Fisher  et  al  [10]  have  shown  experimentally  that  maximum  allowable  value  of  s 
(contractivity  factor)  can  be  1.5  (>  1).  Also  they  have  shown  that  this  maximum 
value  of  s,  for  a  particular  image,  yields  minimum  distortion  between  the  original 
image  and  the  attractor  evolved  through  the  iterative  process  of  the  eventually 
contractive  transformations. 

In  the  next  section  we  have  described  the  mathematical  formulation  of  attractor 
of  PIFS,  where  the  meaning  of  eventually  contractive  maps  has  been  defined. 
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5  Mathematical  Formulation  of  PIFS 


In  this  section  we  have  proposed  a  mathematical  formulation  of  PIFS.  To  make  it 
convenient  we  have  divided  our  tasks  into  three  stages.  Firstly,  it  has  been  shown 
that  the  PIFS  code  {J^)  possesses  a  fixed  point  or  attractor  in  iterative  sequence. 
Secondly,  the  eventual  contractivity  of  the  maps  in  PIFS  setup  has  been  proved. 
Finally,  it  has  been  shown  that  the  given  image  and  the  attractor  are  very  close  to 
each  other  in  the  sense  of  a  chosen  distortion  measure  which  is  root  mean  square 
(RMSE). 

Let,  /  be  a  given  image  having  size  w  X  w  and  the  range  of  gray  level  values  be 
[0,5f].  For  this  given  image  we  can  construct  a  vector  x  whose  elements  are  the 
pixel  values  of  the  given  image  /.  Note  that  there  are  pixel  values  of  L  Thus, 

X  =  (.Ti,  0:2,  ^3,  ...  ,X^2)‘ 

is  the  given  image  where  xi  is  the  pixel  value  corresponding  to  the  (1,  l)th  position 
of  L  Likewise,  let  Xr  be  the  pixel  value  corresponding  to  the  {i^j)th  position  of 
/,  where,  r  =  (i  —  1)  lu  +  j,  1  <  i,  j  <  w. 

In  this  setup  PIFS  can  be  viewed  as  following.' There  exists  an  afRne  (linear),  not 
necessarily  strictly  contractive,  map  for  each  element  of  x  and  this  map  is  called 
forward  map  of  the  element.  In  the  process  of  iteration,  the  input  to  a  forward  map 
will  be  any  one  of  the  elements  of  x  and  the  map  is  called  backward  map  for 
this  input  element.  Thus  for  each  element  of  x  there  exists  a  forward  map  and  an 
element  of  x  can  have  one  or  more  or  no  backward  map(s).  The  set  .F,  of  forward 
maps,  is  called  the  PIFS  codes  of  /. 

Now  let  us  consider  the  set  S  where, 

s  =  {x  I  X  =  (xi,  X2,  X3,  ...  ,X^2)'  ,  0  <  X,-  <  5  }. 

S  is  the  set  of  all  possible  images.  The  given  image  I  is  surely  an  element  of  5 
i.e.  I  G  S'.  The  PIFS  codes  J-  can  be  looked  upon  as 

T  ■.  S  ^  S  . 

The  attractor  of  a  (say),  if  exists  will  also  belongs  to  S.  So,  the  first  task  is  to 
show  the  existence  of  a. 
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Let  /i  be  the  forward  map  for  a  particular  element  ,  where  ri  =  {ii  —  l)w+ji  . 

Also  let  this  element  be  mapped  from  the  element  (  ^2  =  (*2  1)  ^  +  i2  )  • 

Thus  fi  is  the  backward  map  for  xr2  .  Again  is  being  mapped  from 

(ra  =  (ia  ~  l)u;  +  ja)  with  a  forward  map  /2  .  Thus  we  have  a  sequence  of 
maps  for  the  element  x^  as  following. 

^  (*2,^2)  ^  (*3,  is)  ^  •••  (*m,  im);  m  <  {w'^  -  1)  . 

(5) 

The  above  sequence  will  be  stopped  at  (  jm  )  if 

(  im+i,  jm+i  )  =  (  4,  ifc  );  for  A:  =  0  or  1  or  2  or  •••  or  m.  (6) 

The  stopping  phenomenon  of  this  sequence  is  mandatory  as  there  are  finite  number 
(  )  of  elements  in  x  .  Moreover  all  the  elements  of  x  possess  same  type  of 

sequence  in  PIFS  codes.  Thus  it  is  enough  to  show  that  the  element  has  got 
a  fixed  point  in  the  process  of  iteration  and  this  will  lead  to  prove  the  existence  of 
a  (attractor  of  x). 

It  is  clear  from  the  sequence  (5)  that  during  the  iterative  process  the  element 
will  have  a  fixed  point  once  the  element  Xr2  is  fixed.  Again  the  convergence  (to  a 
fixed  point)  of  the  element  X7.3  confirms  the  convergence  of  the  element  Xr2  and 
likewise  for  the  rest  of  the  elements.  Thus  convergence  of  the  last  element  of  the 
sequence  implies  the  convergence  of  the  rest  of  the  elements.  The  convergence  of 
the  last  element  of  the  sequence  is  possible  in  four  different  ways  according  to  the 
stopping  condition  (6). 

An  important  point  to  be  noted  in  this  context  is  the  problem  of  discretization. 
To  get  the  decoded  image  in  an  iterative  process  using  PIFS  codes  one  need  to 
discretize  the  output.  This  can  be  done  in  two  ways.  One  is  discretization  of 
the  output  in  each  iteration.  Another  is  discretization  at  the  end  of  the  iterative 
process.  The  iterative  process  is  stopped  whenever  there  is  no  change  in  gray 
values  in  two  successive  iterations.  To  prove  the  convergence  of  the  elements  in 
four  different  ways  we  have  used  the  discretization  of  the  second  type. 

Case  1  :  m  =  1. 

Here  (  )  =  (  *i,  Ji  )  • 

It  implies  that  (  t'l,  ji  )  is  mapped  into  itself  with  a  map  fi  . 
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Here  fi  =  ai  x  +  bi  ;  0  <  x  <  g  and  0  <  ai  <  1. 

Note  that  in  this  case  the  affine  map  /i  should  necessarily  be  a  strictly  contractive 
map  otherwise  the  element  will  not  converge  to  a  fixed  point. 

If  we  start  with  any  value  (0  <  x  <  g  )  of  (  ij,  ji  )  ,  the  element  will  converge 
to  the  fixed  point  • 

Case  2  :  m  >  0  and  k  —  m 

Here  (  im+l,  jm+l  )  =  {  im,  jm)  • 

It  implies  that  {im^jm)  is  mapped  into  itself  with  a  map  fm  =  o,m  ^  +  bm  ; 
Q  <  X  <  g  and  0  <  <  1  .  Thus  (  jm)  will  converge  to  .  Once 

(  jm)  is  fixed  at  ,  the  element  (  im-i,  jm-i)  will  be  fixed  at 


®m— 1  bm 

1 


+  bm-l 


Qm— 1  bjn _ dm  ^m— 1  ~t~  1 

1  ^m 


In  this  case  the  forward  map  is  ^m — i  ^m“i  x  +  bm-^  ;  0  <  X  <  g  .  Again 
(  *m-i>  jm-1  )  is  fixed  implies  convergence  of  (  im-2,  jm-2  )  with  forward  map 
jm— 2  ”  ^m— 2  X  “b  ^m— 2  j  0  ^  X  ^  ^  ,  at 

^m— 2  ^m— 1  bm  ^m— 2  ®m  bm~~\  "b  ^m— 2  ^m— 1  ^m  ^m— 2  "b  ^m— 2 


Proceeding  in  this  way,  the  fixed  point  of  (  ii,  j\  )  is  found  out  to  be 

ai  02  ...  Om-1  4  (<*1  Q2  •  •  •  «m~2  °2  •  •  •  °m~3  ^m-2  +  •  •  •  -»■  <*1  °2  ^3  +  °1  ^2  +  ^l) 

1  —  am 

Note  that  in  this  case  the  affine  map  /„  should  necessarily  be  contractive  in 
strict  sense.  But  the  rest  of  the  maps  need  not  be  strictly  contractive.  The  eventual 

771 

contractivity,  associated  with  the  element  =  (^i  ?  Ji)  7  will  be  =  JJa*-  . 

i=l 

Case  3  :  m  >  0  and  A:  =  1 
Here  (  ^771-1-1 7  Jm+i  )  (  Ji)  * 

It  implies  that  the  starting  and  the  last  element  of  the  sequence  (5)  is  same.  This 
can  be  looked  as  a  complete  loop  for  the  sequence.  This  case  has  been  solved 
stepwise.  First  of  all  the  case  is  solved  for  m  =  2,  and  m  =  3  .  Then  on  the 
basis  of  these  the  fixed  point  for  the  case  of  general  m  is  solved. 


Case  3(a)  :  m  =  2 

Here  we  have  only  two  elements  viz.  (  ii,  ji  )  and  (  22,  ja)  ■  The  element 
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(  *1,  3 1  )  is  being  mapped  from  the  element  (  *2,  h  )  by  the  affine  map 
f\  —  a\  X  +  hi  ‘  Q  <  X  <  g  .  On  the  other  hand  the  element  (  j2  )  is 
being  mapped  from  (  t'l,  ji  )  by  the  affine  map  /j  =  <12  2;  +  62  !  0  <  x  <  g  . 

(  *1)  ji  )  ^  (  *2!  32)  ^  (  *1,  ii  )  • 

Let  x  be  the  starting  value  of  (  ii,  ji  )  and  y  be  the  starting  value  of  (  22,  72  )  • 
After  first  iteration  the  values  of  (  ii,  ji  )  and  (  12,  j2  )  will  be  Oi  y  +  hi  and 
02  X  +  62  respectively.  Again  after  second  iteration  these  will  be  oj  02  x  +  oi  62  +  hi 
and  02  ai  y  +  02  ^1  +  ^*2  respectively.  Proceeding  this  way  after  infinite  (practi¬ 
cally  large  but  finite)  number  of  iteration,  the  fixed  point  of  (tj  ,  ji)  and  (22  , 
will  be  independent  of  x  and  y.  The  Coefficients  of  x  and  y  after  N  (even)  itera¬ 
tions  will  be  (ci  02)^^^  which  tends  to  zero  as  N  tends  to  infinity.  The  fixed  points 
of  (  *1,  ji  )  thus  will  be 

Qi  1>2  d" 

1  -  oi  02' 

The  same  for  the  element  (  22,  72  )  will  be 

02  61-1-62 

1  —  Oi  02  ’ 

Note  that  both  the  maps  need  not  be  contractive.  Moreover  the  eventual  con- 
tractivity  associated  with  the  element  x^  is  (oi  02)  which  should  be  less  than 
one. 

Case  3(b)  ;  m  =  3 

Here  we  have  three  elements  viz.  (  z’l,  ii  )  ,  (  12,  j2  )  and  (  23,  jg)  .  These  three 
elements  are  making  a  complete  loop  in  the  sequence.  The  sequence  of  forward 
and  backward  maps  is  as  follows; 

(  *1)  3l)  ^  (  *2,  i2  )  ^  (  *3,  is  )  ^  (  h,  3l  )• 

Taking  the  starting  values  of  three  elements  as  x,  y  and  z  and  proceeding  as  case 
S(o)  we  have  the  following  results. 

The  fixed  point  of  (  ii,  ii  )  will  be 

Ol  (  02  63  +  62  )  -f  61 

1  —  Ol  02  O3 
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The  element  {  2*2,  32  )  will  converge  to 

^2  (  ^3  ^1  +  ^3  )  +  ^2 

1  —  fli  a2  03 

The  fixed  point  of  (  «3>  )  will  be 

4*3  (  62  +  61  )  +  63 

1  —  Oi  Ct2  dz 

Here  also  the  maps  need  not  be  contractive  in  the  strict  sense.  The  eventual 
contractivity  will  be  (ai  a2  a^)  in  this  case. 

Case  3(c)  :  General  m 

Here  we  have  m  elements  which  are  making  a  complete  loop  of  sequence.  It  is  clear 
from  case  3(a)  and  case  3(b)  that  all  the  elements  of  this  sequence  will  have  a  fixed 
point  after  a  large  but  finite  number  of  iteration.  Also  the  affine  maps  which  are 

used,  need  not  to  be  contractive.  In  particular,  in  this  case  the  element  (  Zi,  ji  ) 

will  converge  to 

ai  (a2  {  •••  (Qm-2  (^m~l  bm  +  ^m-l)  +  bm-2)  +  »»»)  +  62)  +  61 

1  a\  a2  •  •  •  dfii 

m 

Also  the  eventual  contractivity  for  the  element  is  5^^  =  ][[  • 

i=i 

Case  4  :  m  >  0  and  0  <  k  <  m 

Here  (  jm+\  )  =  (  4,  jk)  ,  where  k  =  2  or  3  or  ...  or  m  -  2  . 

Without  loss  of  generality  say,  l<fc  =  mo<m  —  1. 

This  case  can  be  viewed  as  mixture  two  cases.  Taking  (  jmo  )  the  starting 
element,  a  complete  loop  of  sequence  can  be  formed  with  rest  of  the  elements. 
Thus,  one  can  find  the  fixed  point  of  this  element  as  it  is  nothing  but  case  3.  Once 
the  element  (  j^o  )  is  fixed  then  the  fixed  point  of  the  original  starting 
element  (  ii,  ji  )  can  be  found  out  by  using  ca^e  2,  Like  all  the  previous  cases 

m 

the  eventual  contractivity,  in  this  case,  will  be  Sr^  =  a*  . 

t=i 

The  above  stated  four  cases  provide  the  fixed  point  of  the  PIFS  codes  T.  Thus  for 
a  very  large  positive  number  N ^  we  have  a  ^  o  ^  5  .  So, 

lim  :F^(£)  =  a,  V  o  e  5  (7) 

jV—foo  ^  ' 

Also, 
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^{a)  =  Jim  V  o  e  5. 

=  lim  :r^+^(o);  Woes. 


Note  that  for  each  element  there  will  be  a  sequence  of  the  form  (5).  This  sequence 
will  follow  any  one  of  the  above  mentioned  four  cases.  Thus  for  each  element  there 
will  be  a  sequence  of  forward  maps.  The  contractivity  factor  associated  with  this 
element  will  be  the  product  of  all  the  scaled  parameters  (  a,-  )  of  the  forward  maps. 

The  next  task  is  to  define,  mathematically,  the  eventual  contractivity  of  the  PIFS 
codes.  In  this  context  we  are  stating  the  following  theorem. 

Theorem  1  :  Let  ^  be  the  PIFS  codes  and  S  be  the  set  of  all  possible  images. 
For  every  x  and  y,  x^y^eS3N>0  and  0  <  s  <  1  such  that 

I  I  <  s  I  X  -  j/  1  ;  W  p,  q  >  N  . 

Proof ;  Using  equation  (7)  we  have  for  a  very  small  positive  number  e  >  0  ,  3 
a  large  positive  number  such  that 

p  >  Ni  =>  1  -  a  I  <  ^  ;  V  X  e  5. 

Also,  3  another  large  positive  number  N2  such  that 

q  >  N2  ^  \:F’>{y)  -  a\  <  i  ;  W  y  £  S. 

Thus  for  X  £  6  S  , 

1  ^'’(21)  -  I  =  I  -  a  +  a  -  jf^(£)  |. 

<  1  ^’’(^)  -  o|  +  I -«l- 

<  e  ;  where  ,  p,  q  >  N  =  Max  {Ni,  N'i). 

Thus,  for  X  7^  y,  3  a  large  number  N  and  0  <  s  <  1  such  that 

p,  q  >  N  I  J^{x)  -  J^’^iy)  I  <  5  I  X  -  y\.  Q.E.D. 

Once  the  eventual  contractivity  of  the  PIFS  codes  ^  has  been  proved  the  last  task 
is  to  show  that  the  given  image  and  the  fixed  point  of  are  very  close  to  each 
other. 
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Theorem  2  :  Let  a  be  the  fixed  point  of  the  PIFS  codes  T  and  x  be  the  given 
image.  Also  let  “d"  be  the  given  distortion  measure.  Under  this  setup,  if 

d(  X  ,  J~{x))  <  a 

then 

d{  X  ,  a)  <  - - - . 

1  ~  ^max 

Where  Smax  =  Max  {si,  S2,  ...  ,5u,2}  ,  s,-  being  the  eventual  contractivity  of 
the  ith  element  of  x  . 

Proof :  Let 

—  ~  (^l>  ^2)  ••• 

:f{x)  =  (xj,  ...  ,x;;3)' 

y^ix)  =  (g,  g,  ...  ,f^)' 

=  (xl,  ... 


In  the  PIFS  scheme,  the  given  image  is  partitioned  into  range  blocks  TZj  of  sizes 
6x6.  So,  there  are  n  =  ( range  blocks  each  having  6^  pixel  values.  Pixel  values 
are  nothing  but  the  elements  of  x.  Also  the  distortion  measure  “d”  is  Root  Mean 
Square  Error  (RMSE).  “d”  ,  defined  on  5,  is  as  follows 


\ 


1=1 


V  e  S. 


Where,  d  (  u,-  ,  u;  )  =  |t/,-  -  . 

Now, 
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X  ,  T{x)  )  _  ^  c?  (x;  ,  fi{xj))\ 

[a:,-  is  being  mapped  from  Xj  with  forward  map 

\  ^  (^*  5  ^t)* 

'  t=l _ _ 

w'^ 

-  \  Maxi\xi  - 


^  H  (<^)^  vvhere,  a  =  Max 


iS{l f2,...,w^y  \^i  ^i\‘ 


Again, 


^  ^  \  ^  ^  )  /i|t(a'’i|;))  ,  [as  =  n6^] 

\  t-1  j=l 

where,  .r,,,-  is  the  jth  pixel  value  of  ith  range  block  and  is  the  forward  map 

associated  with  x,,,.  which  is  being  mapped  from  x,,,  ,  the  element  of  range 
block.  ^ 


1  ^  1  62  ^ 
\  S  ^  ^  ^  ^  ,  x~j}). 

It  implies  that 


n  ^  ^  ^  ^il'  ’  ^  O! 

«=i  "  ;=i 


Again, 


)  =  i  t  ^  I;  ,(  ^ 

\  t=l  j=l 

Note  that  the  size  of  the  range  block  and  the  contracted  domain  block  (from  where 
t  IS  range  block  is  being  mapped]  is  same  (Section  3).  Moreover  the  number  of 
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range  blocks  and  the  number  of  matched  contracted  domain  blocks  is  same  as  there 
is  only  one  matched  contracted  domain  block  for  each  range  block.  Also  for  each 
element  there  is  an  eventual  contract! vity  factor  and  Smax  is  the  maximum  of 
these  factors. 


d{J='{x)  ,  J^{x))  <  Sr, 


1  A  1  ^ 


12  ^2 


”  /=i  ^  i-=i 


<  Smoi  o  •  [By  (8)] 


Similarly, 


,  3^(1)  )  =  i  E  n  i:  <i(  ,  /,|,(i5))  • 

I  1  A  1  ^  ^ 


^  ^max  \  ^  12  ^  d{  XkH  ,  . 


„  ^  A2 
"  ^=1  ^  A;=l 

1  A  1  ^ 


A  ■»  1  - 

=  Smax  \  ;r  zJ  ^  A  /fciK^^pk)  >  A-|i(^pk))  • 

/=i  ®  k=i 


<  si 


[where  Xk\i  is  being  mapped  from  Xp\g 
with  forward  map  fk\t  ] 

n  2  6^ 

^  ^  ^  ^p\i  ’  ^pk)  • 


<  a  .  [By  (8)] 

So,  finally  we  have 

d  {  X,  a)  =  d  {  X,  )  ;  'i  o  ^  S  . 

=  d  {  X  ,  {x)  )  ;  [by  Theorem  1] 

<  d{x  ,  J^{x))  +  d{J^{x)  ,  r^{x))  +  d{J^\x)  ,  .F"(£))  +  .... 

Q  +  Smax  Ol  +  Q;  +  .  .  . 

CC  (  1  +  Smax  +  'Smai  "b  •  •  •  )  • 

=  •  Q-  E.  D. 

In  the  next  section  we  have  presented  the  experimental  results  in  support  of  the 
mathematical  formulation  of  PIFS. 
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6  Experimental  Results 


In  the  context  of  PIFS,  a  technique  for  fractal  image  compression  using  Genetic 
Algorithm  (GA)  has  been  proposed  by  Mitra  et  al  [11,  1].  Using  this  technique,  the 
PIFS  codes  for  256  x  256,  8  bits/pixel  “Lena”  image  and  “LFA”(Low  Flying  Air¬ 
craft)  image  has  been  found  out.  Both  images  have  been  reconstructed  iteratively 
starting  from  different  images.  In  particular  we  have  used  a  “Blank”  image  (having 
all  the  gray  values  zero),  “Rose”  image,  “Lena”  image  and  “LFA”  image  as  start¬ 
ing  images.  At  the  time  of  reconstruction,  RMSE  (distortion  measure)  between 
two  successive  iterations  has  been  computed.  In  all  the  cases,  computed  values  of 
RMSE  are  gradually  decreasing  with  the  increasing  iteration  number.  It  has  been 
found  that  the  PIFS  codes  almost  achieved  a  fixed  point  after  ten  (10)  iterations 
in  both  the  images.  Finally,  the  distance  (RMSE),  as  it  is  expected,  between  the 
given  image  and  the  attractor  is  found  to  be  very  small.  The  distances  are  found 
to  be  on  an  average  7.75  and  11.44  for  “Lena”image  and  “LFA”  image  respectively. 

The  computed  values  of  RMSE  have  been  presented  in  tabulated  form  in  Table  1 
and  Table  2.  The  notation  Aij  used  in  the  tables  implies  the  RMSE  value  between 
ith  and  jth  iterations.  Thus  we  have 

A,-,-  =  d  (  ,  P{o)  )  , 

where  JF  is  the  PIFS  codes  used  and  o  is  the  starting  image.  As  we  have  stopped 
the  process  of  iterations  after  10  iteration,  so  the  value  of  Ao,io  provides  the  dis¬ 
tance  between  the  attractor  and  the  given  image.  The  value  of  Aoi  will  provide 
an  approximate  value  of  a  if  the  starting  image  is  the  given  image  itself.  The 
approximate  values  of  a  are  found  to  be  7.18  and  10.67  for  Lena”  and  LFA” 
images  respectively. 


Table  1:  RMSE  values  between  successive  iterations  using  PIFS  codes  of  “Lena” 


image 


Starting 

image 

^0.10 

^0.1 

^1.2 

•^2.3 

^3.4 

^4.5 

^5.6 

^6.7 

^7.8 

^8.9 

^9.10 

Lena 

7.73 

7.18 

2.91 

1.12 

0.45 

0.24 

0.15 

0.10 

0.07 

0.04 

0.03 

Blank 

7.86 

X 

43.93 

32.05 

19.72 

11.30 

6.25 

3.48 

1.97 

1.13 

0.66 

Ro^e 

7.74 

99.96 

62.42 

37.63 

19.88 

9.45 

4.13 

1.91 

0.92 

0.48 

0.28 

LFA 

7.74 

66.47 

41.59 

23.07 

11.26 

4.97 

2.13 

0.97 

0.48 

0.27 

0.17 
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Table  2:  RMSE  values  between  successive  iterations  using  PIFS  codes  of  “LFA” 


image 


Starting 

image 

^0.10 

^0.1 

^1.2 

^2.3 

^3.4 

^4.5 

-^5.6 

^6.7 

^7.8 

^8.9 

-^9.10 

LFA 

11.42 

10.67 

4.32 

1.56 

0.63 

0.34 

0.23 

0.16 

0.12 

0.09 

0.06 

Blank 

11.48 

X 

53.78 

36.51 

21.87 

12.70 

7.39 

4.28 

2.46 

1.42 

0.84 

Rose 

11.45 

143.09 

102.65 

43.39 

22.34 

10.47 

5.18 

2.74 

1.53 

0.89 

0.56 

Lena 

11.43 

64.95 

45.50 

27.23 

13.38 

5.66 

2.36 

1.09 

0.60 

0.38 

0.26 

To  judge  the  validity  of  the  PIFS  technique,  the  original  and  reconstructed  images 
have  also  been  checked  visually.  Figure  2,  Figure  3  and  Figure  4  are  showing  the 
original  images  of  “Lena”,  “LFA”  and  “Rose”  respectively.  Figure  5  and  Figure  6 
are  showing  the  reconstructed  images  of  “Lena”  and  “LFA”  respectively.  In  both 
the  cases  the  starting  image  is  the  “Rose”  image.  Both  the  reconstructed  image  are 
found  to  be  of  good  quality.  Also  the  compression  ratios  are  found  to  be  10.5  and 
5.5  for  “Lena”  and  “LFA”  images  respectively,  using  GA  based  PIFS  technique  [1]. 
Note  that  the  compression  ratio  is  other  factors  like  image  size,  range  block  size, 
bits  per  pixel  and  the  gray  level  variation  present  in  the  given  image. 

The  next  section  concludes  the  present  article. 


7  Conclusions 

The  upper  bound  of  the  difference,  between  the  original  image  and  the  attractor 
evolved  through  its  PIFS  codes,  is  almost  same  as  that  in  the  IFS  set  up  (Collage 
theorem  [2]).  The  only  difference  is  the  contractivity  factor  which  is  eventually 
contractive  in  the  case  of  PIFS  compared  to  the  strict  contractivity  in  IFS.  The 
present  article  provides  an  elaborate  and  direct  proof  for  the  existence  of  the  at¬ 
tractor  and  the  closeness  of  the  attractor  to  the  given  image  in  the  partitioned  IFS 
scheme. 

In  PIFS  technique  the  estimates  of  all  the  range  blocks  are  obtained  assuming  the 
self  similarities  present  in  the  given  image.  The  domain  block  which  is  most  similar 
to  a  range  block  is  named  as  appropriately  matched  domain  block  for  that  range 
block.  The  similarity  between  the  range  block  and  the  domain  block  is  measured 
by  RMSE  (or  MSE).  Thus  the  efficiency  of  PIFS  technique  depends  on  two  facto^?. 
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The  first  one  is  the  eflficiency  of  the  distortion  measure  and  second  one  is  the  extent 
of  similarity  present  in  the  given  image. 

RMSE  being  a  global  measure  has  its  own  limitations.  In  this  context  one  can 
think  of  a  better  and  reliable  measure  which  can  make  the  PIFS  technique  more 
efficient.  Regarding  the  second  factor,  it  may  so  happen  that  there  is  hardly  any 
domain  block  which  is  appropriately  matched  with  the  concerned  range  block.  In 
other  words,  the  domain  block  most  close  to  a  range  block  in  the  sense  of  similarity, 
may  provide  a  quantitatively  large  distortion  measure.  This  may  lead  to  inefficient 
coding.  This  problem  can  be  viewed  as  a  limitation  of  the  PIFS  image  compression 
technique.  Thus,  there  is  enormous  scope  for  suggesting  the  specific  theory  which 
can  make  domain  blocks  more  close  to  the  range  block  in  the  sense  of  similarity 
and  also  with  less  distortion. 

The  theory  of  PIFS  technique  is  also  applicable  to  one  dimensional  signals.  In  This 
context  the  technique  has  already  been  applied  to  code  the  fractal  curves  [7]  and 
EEG  signals  [13].  Computational  time  and  hence  the  computational  cost  of  the 
PIFS  scheme  for  gray  level  image  compression  is  quite  large.  So,  researches  are 
carried  out  to  reduce  the  computational  time  as  well  as  computational  cost  [1,  14]. 
Also  the  PIFS  codes  may  be  used  in  other  aspects  of  image  processing. 
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Table  Captions 


Table  1.  ;  RMSE  values  between  successive  iterations  using  PIFS  codes  of 

“LENA”  image  image. 


Table  2.  :  RMSE  values  between  successive  iterations  using  PIFS  codes  of  “LEA 
image  image. 


Figure  Captions 


Mappings  for  IFS  scheme 

Mappings  for  domain  blocks  to  range  blocks  in  PIFS  scheme 
Original  “Lena”  image 
Original  “LFA”  image 
“Rose  image 

Decoded  “Lena”  image  after  10  iterations  using  the  starting  image  as  Figure  4 
Decoded  “LFA”  image  after  10  iterations  using  the  starting  image  as  Figure  4 


